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An essential feature of ultra cold atomic gases, which are intensively studied today, is there finite volume. From 
the point of view of the well known Lee - Yang theorem the phase transitions in the finite size systems are impossible 
0. An influence of finite sizes on the coherent properties of a dilute Bose gas in a trap is analyzed in the present 
work. The trap is described by the stationary external potential U { r ) which confines the system in a finite volume. 
The external potential is supposed sufficiently smooth, so that the chemical potential fi = ng \s much larger than the 
distance between one-particle levels of the external potential U (1^), n is the average density of the gas and g is the 
coupling constant characterizing the repulsive interaction between particles. Here and below we choose the system of 
units in which the particle mass and Planck constant are equal to unit m = 1, h = 1. 

The problem of the fluctuations of the condensate in the finite system with the fixed total number of particles 
due to the long- wave phonon excitations has been considered in the work 2]. The anomalous dependence of these 
fluctuations on the volume of the system was obtained there. The phase fluctuations of the condensate in a highly 
elongated traps have been considered in series of works, for example 0. The mechanism of fluctuations of the 
condensate in the flnite Bose systems connected with the fluctuations of the global phase has been considered in a 
series of works, for example, 0, 0i 0- these works the correlation function K\\ =< ip (~r,t) ip^ {~r',t') >, ip and 

'ip'^ being the annihilation and creation operators of Bose particles, has been calculated and it has been shown that 
this correlator decays exponentially for t =\ t — t' \ — > oo due to the phase diffusion. 

In the present work we consider the last mechanism of large fluctuations in the flnite Bose systems connected 
with the fluctuations of the global phase and the existence of the zero mode, i.e., the zero energy excitation mode, 
corresponding to the invariance of the action related to a change of the global phase j^. The method which we use 
differs from the methods of the known works Q and gives the possibility to consider quantum fluctuations 

and dissipation in the unifled technique. 

Due to existence of zero mode separated by the gap from the other excitation spectrum, the phase fluctuations 
in flnite systems have an infrared frequency divergence independent of the space dimension. On the other hand, in 
the limit of inflnite sizes this infrared divergency disappears. As a result, the properties of the correlation functions 
in the finite and infinite systems are essentially different. Thus, the calculations of normal and anomalous 

K± =< ip (~r*, t) "0 (~?^', t') > correlation functions at low temperatures show that correlator K\\ decays exponentially 
in time for r =| t ^ t' \ — > c» and the anomalous correlator K±^ vanishes for all r. The reason of such behavior 
of these correlators is the infrared frequency divergency of phase fiuctuations in the finite systems. At that, due to 
the disappearance of this divergency in the case of the 3D system of infinite sizes L — > cxd, <N>— > cxd, where L is 
the characteristic size of the system and <N> is the total number of particles, the normal and anomalous correlators 
restore their usual properties, i.e., these correlators are equal to no where no is the condensate density. Note that the 
anomalous averages are considered in the sense of quasi-averaging, i.e., the infinitesimal term j^if) -\- breaking the 
global gauge invariance is added to the Hamiltonian and the source j+ tends to zero at the end of the calculations. 
This infinitesimal term results in the existence of the infinitesimal gap in the energy of zero mode. The infinitesimal 
term will not be written evidently but the existence of it will be implied always. The average < ■0 ( r , t) > , which 
usually considered as the order parameter, is calculated too and being considered in the sense of quasi-averaging gives 
the zero value. Note that all these averages are calculated in the supposition of the thermodynamical equilibrium and 
are accomplished over the equilibrium density matrix. 
The action of a nonideal Bose gas in a stationary trap is 



V'} (1) 

The contour of the integration over the time is the two-time contour which was introduced in the work . We use 
the Keldysh - Schwinger technique Q , for the consideration of the system at low but finite temperature and in the 
real time. 

The anomalous and normal correlation functions can be written in the form of the functional integral as 



S = idt 



idt 
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-U{~r ) 



Kj_^\\ (Vi,ti;T*2,t2) = J pDpDif exp {iS} p (xi) p {x2) exp {i {'f {xi) ± (f{x2))} 

where we denote x = {~r^,t), the complex field ip is written in the modulus-phase representation -0 = pe^^, here 
and further the sings and (-) relates to the anomalous and normal correlators respectively. The under integration 
expression for the anomalous and normal correlators can be represented as an exponent with the actions Sj(±-), Sjuo 
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respectively 



K±^\\{l^i,ti;1^2,t2) = j DipDnp{xi) p{x2)exp{iSj(±,ii)} (2) 

where Sj{±,\\) are 



5^(^,1) = fdtj d''r|j(^'ll)(^ + n^-in(V(^)'-i(V/9)'-|n2 + n(/x-?7(y))| 
Here we denote n = p^, = dtf and 

It is convenient to introduce new variables p ^ p and ip ^ pip that gives 

S = jdtj d^r|2pC-^%-^(Vp)'-|n2 + n(/x-C/(r))| 



(3) 



(4) 



where 



and A = ^. 



p 



Denote the orthonormal system of eigcnfunctions and eigenvalues of the Hamiltonian as rjx and £| (A) , corre- 
spondingly, A is a quantum number of the eigenfunction r]\, and 

HiVx = £i (A) V\ 

The spectrum of the operator ffj has the following properties: firstly, it is a discontinuous one and A can be 
considered as an integer, secondly, the eigenvalues of this operator have the zero or positive value (A) ^ 0, thirdly, 
the zero value (0) = exists for any configuration of the field p (T*, t) and it is separated from the other spectrum 
£5 (A) (A 0) by the gap which goes to zero when L — > 00. 

The system of eigcnfunctions rjx can be chosen as an orthonormal system 

< V\' I V\ >= 

A wave function of the zero mode £j (0) = has the form 

r}o{'r,t) = -^p{'r,t) (6) 



The existence of the zero mode can be seen directly from the action of the operator (5) on the wave function (6). 
The eigenfunction expansion of ^ is 

i{'^,t) = Y,nx{'f,t)ax {t) 

A 

and the zero mode in this expansion corresponds to the value A = 0. Below we separate the zero mode in the 
eigenfunction expansion 

^ (r , i) = ^ r?A t) ax {t) = vo t) ao {t) + ^' , t) (7) 
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where 



e{r,t) = Y,Vx{^,t)ax{t) (8) 



As a result we obtain 



dt<-^ao (i) 



Hi / 
-n +n(/i- 



(9) 



Note, that the first term in the action S (9), as a result of the integration over the zero mode {t) in the functional 
integral for correlators (2), gives the conservation law of total number of particles , namely, the functional 5-function 

. This ^-function gives the constraint for the density fluctuations, and this constraint has to be taken into 

account when the subsequent integration over the field p is accomplished after the integration over ao [t). Thereupon, 
the simplest way to take into account the constraint connected with the conservation law of total number of particles 
is the integration over ao {t) after the integration over the density fluctuations, i.e., over the fleld p. 

For the integration over the field p the saddle-point approximation can be used. In this approximation the value of 
the field p should have the most symmetrical form and, thus, does not depend on the time t. The field p is represented 
as p (7^, t) = po (r') + Sp (V. t), where po C^) is the saddle-point configuration and 5p {T^, t) is the fiuctuation. The 
saddle-point configuration po (V) is defined by the Eq. 



-^V^'P + Uir)+gpl po = 
The variation of S can be considered up to the second order of fluctuations 5p and can be represented in the form 



6S = j^dtlao {t) Wj^dtjd'^r !^26p^' - ^'H^'^^^' - SpH^^dp 



(10) 



where 



i^f -V^ + 



Wo 

PO 



and 



Hf) = -\v^-Ui + Zgpl + U{r) 

The value 5N is 5N {t) = J d''r2po ( r") Sp {'r, t) and iVo = / d'^rpl {'r). 
The integration over 5p in the Eqs. (2) for the correlators K±_ and K\\ gives 

-f_L,|| (^i,^2;ii,i2) =Po(^i)po (^2) j r'Cexpj^S'j [^]} 

where 
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As it will be seen later the main contribution to the functional integral over the field qq (t) gives the low frequency 
configurations, such that their frequencies are much smaller than the frequencies of the field . Due to this fact 
the terms proportional to oq^' can be omitted in the action (11). Moreover, taking into account only the slow 

configurations of the field ^' which have the wave vectors smaller than y//I we can replace the term (^HjP^^ by 

(~r^ — 1^') and, thereby, rewrite the action (11) as 

^S., Kl = / (ff - f Hf 'f + J + } ,12) 

The integration over uq (t) and (^,i) gives the following result for the correlators 

-^±,11 (^1, ^2; ti,t2) = po ( r*i) po (^2) exp {i5'j_^|| } (13) 



-1 fdu; 



D 



?^ (^) + Dg^ (iu) ± (c) + D^^^ (c.)) cos (ujt) 



TT {D'li {^1,^2-,^) + D22 (^1,^2;^^) ± [D[2 (^l,^2;w) +i>21 (^1,^2;W)]C0S {lJt)} 
ZTT 

where £'^''^(0;) is the Keldysh-Schwinger Green function in "four-angle" representation ;3j for the field ao (i) and 
D'^ p{r I, r 2; w) is the Keldysh-Schwinger Green function for the field ^' {r ,t), indices a, (3 = 1,2 means the upper 
and the lower branches of the Keldysh-Schwinger time contour. As before, the upper and the lower signs relate to 
K± and ify respectively. Eq. (14) can be rewritten in the form 

^±.11 ^ / (u;) (1 ± cos icjT)) 2 / (A; cj) r^, (V 1) (^2) (1 ± cos {u;r)) (14) 

The Eq. (14) is obtained from the Eq. for S'j_j| (13) using the following property of Keldysh-Schwinger Green 
fimctions in "four-angle" representation 

Dii [lj) + D22 {uj) - D12 (uj) - D21 (iu) - 

and the following connections between non-diagonal components of these Green functions in " four-angle" representa- 
tion, "kinetic" Green function Dk, as well as the retarded Dji and advanced Da Green functions 

D12 (w) + D21 {oj) = 2Dk (w) = 2coth (^) [Dr (w) - Da (w)] 



The last Eqs. take place in equilibrium states. If the relaxation processes are not taken into account the Green 

)(0) 

'k 



functions Dr,a (f^) and D^^ {uo) have the form 



^^^^(^^-)- .^-2,e'(A)±.^ ^ i.^")(.) = 2..coth(|,)..5n(.).4-^] (15) 

The regularization of the first integral in the Eq. (14) by means of the change 5 [o;^] — > 5 [cj^ — £q] and then the 
transition to the limit eo ^ results in 



2 T^mA „ . 2 2 

1 TT — 9- ; '-'II = «TT--' 

eo^o elJ' " iVo 

, at that, the first integral in the Eq. (14) gives the main contribution to •S'j.,!!, thus, it gives the infinite contribution 
to S_L^ and the second integral gives the finite contribution to S±_ due to the existence of the gap in the spectrum of 
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the 5' excitations. Obviously, that these Eqs. result in the zero value for the anomalous correlator K±^ (t) = and 
the exponential decrease of the normal correlator K\\ (r) — Pq exp (^—^Tpr^^ for the sufhciently large r. 
Note, that in the limit of the infinite size of the system L— > oo the values S_\_, Su transfer to 



^ _ j 2 /■ duid'^k 



Dk (uj,~k^ {I ± cos (ujt)) (16) 



where 



Dk (^i^, ^ ) ~ 27ri coth sign (lu) pS [w^ 



2 ;.21 



C^k 



In this case due to the integration over the momentum k there is no divergency of the integral over the small 
frequencies and with an accuracy to the small gas parameter in the dilute gas the anomalous and normal correlators 
takes the usual form = K± = in 3D case. This fact can be obtained easily by the calculation of the integral 
(16), more over, this calculation gives the well known vanish of K± in 2D case for a nonzero temperature and in ID 
case for any temperature. 

The singularity of the expression for d'^^P {uj) strongly changes by the relaxation of the phase fluctuations. To obtain 
the action taking into account the relaxation, the terms proportional to the fourth power of Sp, i.e., g{Sp)'^ should 
be considered. Integrating over the rapid fields Sp and i.e., the fields having the momentums k > T/c, where c 
is the sound velocity and the temperature T is supposed being much smaller than p, we obtain the additional term 
SpYjSp in the effective action for the slow fields. In the case of the dilute Bose gas the self-energy part S is determined 
by the diagram of the second order in the coupling constant g and has three internal particle lines corresponding to 
Green functions of the rapid 5p field Gsp,Sp — —i < Tt5p{x)5p{x') >, where Tt means t-ordering of operators. The 
additional self-energy part for the slow 5p field can be written as S {xi,X2) = g^ {Gsp,5p {xi,X2)) and, namely, this 
term describes the relaxation processes. In the case of the dilute Bose gas only the imaginary part of this self-energy 
part can be taken into account, the real part can be omitted due to its smallness compared with the chemical potential. 
The self-energy part S for small frequencies w — > is calculated as ^ fl^^'^'^ ^^"^ ^^'^ kinetic part T,k [1| 

takes the form ^ . Note that the internal momentums which give the main contribution to the diagram for 
S are pi P2 ~ T / c. After the integration over 8p and taking into account the additional term SpUSp the function 
r (w) is calculated as T {lu) - -^w^T^. 

After the transformation to the action being written in terms of slow fields "0 and containing the relaxation term 
the transition to the modulus-phase variables for the slow fields can be made. After that the subsequent scheme of 
the transformations is analogous to those which has been made above, i.e., without consideration of the relaxation 
term. This gives the expressions for in the form (14) where the phase correlators d'^^\ [uj) and ^ (w. A) have 
to be changed by 

In the case of the low temperatures T « p and for small frequencies to —t Q the relaxation V [lo) can be calculated as 

r {ijj) ^ uj ^'^^'i g^ , where c is the sound velocity. Note, that the internal momentums which give the main contribution 
to the diagram for the relaxation T (lu) are pi ^ P2 ^ T/c. At that, the kinetic Green function has the form 

The kinetic Green function {uj) for small frequencies w — > can be written as D^-* {uj) ~ -»A'^r(i^) ^^^^i (^y) ~ 
—i ^6^2 g'^T^ : i-e., it has the nonintegrable singularity for — s- 0, in contrast to the correlator Dj^ (w, A). Thus, the 
integral for S± (14) diverges at small frequencies as I/lo. Due to the divergency of S± the anomalous correlation 
function Kj_ vanishes Kj_ (r) = 0. The normal correlator i^y (r) for sufficiently large r decays exponentially and 

has the form Ki\ (r) = p§ exp ( — ^2p^T ^-^^ I I ) • Note that nonzero value of the relaxation F (w) for small 
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frequencies is obtained only in the case when the deviation of the excitation spectrum from the hnear one for small 
momentums has the positive sign, and when this sign is negative the relaxation T (lu) has zero value. 

In conclusion, we note that the calculation of the order parameter < t/j > is analogous to the calculation of K± and 
gives the lack of the order parameter <?/;>= in the finite Bose systems as well as vanishing of K±. On the other 
hand, in the 3D case in the limit L— > oo the anomalous correlator Kj_ as well as the normal correlator K^^ and the 

order parameter < ijj > recover their usual values K± = K^^ = ^< ^ >^ = Pq. 
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